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BaH (and its isotopomers) is an attractive molecular candidate for laser cooling to ultracold temperatures 
and a potential precursor for the production of ultracold gases of hydrogen and deuterium. The theoretical 
challenge is to simulate the laser cooling cycle as reliably as possible and this paper addresses the generation 
of a highly accurate ab initio potential for such studies. The performance of various basis sets within the 
multi-reference conhguration-interaction (MRCI) approximation with the Davidson correction (MRCI-I-Q) is 
tested and taken to the Complete Basis Set (CBS) limit. It is shown that the calculated molecular constants 
using a 46 electron Effective Core-Potential (ECP) and even-tempered augmented polarized core-valence basis 
sets (aug-pCVnZ-PP, n = 4 and 5) but only including three active electrons in the MRCI calculation are 
in excellent agreement with the available experimental values. The predicted dissociation energy Dg for the 
state (extrapolated to the CBS limit) is 16895.12 cm“^ (2.094 eV), which agrees within 0.1% of a 
revised experimental value of <16910.6 cm“^, while the calculated re is within 0.03 pm of the experimental 
result. 


I. INTRODUCTION 

Recently, it has been proposecP that ultracold hydrogen 
atoms may be formed by the photo-dissociation of laser 
cooled hydrides, provided the energy of the fragmenta¬ 
tion products is much smaller than the average thermal 
energy of the parents. At the limit of zero kinetic en¬ 
ergy release, the velocities of the hydrogen atoms will 
match that of the parents, forming a Maxwellian veloc¬ 
ity distribution corresponding to a much lower hydro¬ 
gen atom temperature than the original molecular gas 
(T// = ^-^" ). Quantum chemical calculations of the 
transition dipoles and Franck-Condon (EC) factors, using 
the post Hartree-Fock methods Complete Active Space 
Self Consistent Field (CASSCFjPand MRCP, confirmed 
that BaH was a very good candidate for demonstrating 
this kinematic effect as the barium atom has a consid¬ 
erable mass and the molecular radical is amenable to 
laser cooling. In the earlier study a very small (triple 
zeta quality) two electron basis set for barium was used 
with the inner 52 electrons replaced by an ECP. To im¬ 
prove agreement with the experimental energy levels of 
the atom, an ^independent Core-Polarization Potential 
(CPP) was added to the calculation. Although the re¬ 
sulting basis set simulated the reported experimental FC 
factors very well, there were clearly issues with the cal¬ 
culation of the ground state equilibrium bond length 
(experimental valu^ 2.23188651(19) A) and the dissoci¬ 
ation limitPf^. This paper attempts to produce a global, 
high-quality ab initio potential for the ground state 

of BaH. The heavier familial hydrides YbH and RaH are 
also attractive precursors of ultracold hydrogen but, un¬ 
like BaH, there is currently a lack of reliable theoretical or 
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experimental data to confirm their suitability. Hydrogen 
forms the basis of many spectroscopic tests of fundamen¬ 
tal symmetr ieJ^, for the evolution of the Univers^^^ and 
in metrologj®^ (for example, the measurement of the 
Rydberg const anlfSl and the proton radiuJi^. 

BaH is naturally a heavier cousin of BeH, the simplest 
stable neutral open shell diatomic. BeH possesses just 
five electrons and h as recently been the focus of two im¬ 
portant paper^iSIi^ regarding the simulation of its ground 
potential energy curve. The firslfi^was a detailed ab ini¬ 
tio study using very large correlation consistent basis sets 
(aug-cc-pCVnZ) up to seven zeta (n = 7) in character 
(indeed so large that the MOLPRCP^ ab initio package 
could not use the largest 1 = 1 functions in the basis set 
and the effect of excluding them was determined by a 
separate calculation). When determining the static elec¬ 
tron correlation beyond Hartree-Fock level, the CASSCF 
wavefunction included the excited 3s and 3p orbitals on 
the Be atom within the active space despite this being 
a single electronic state calculation. To avoid the size- 
extensivity issues inherent in MRCI (for dynamic elec¬ 
tron correlation) the MR-ACPJ^i^multi-reference aver¬ 
aged coupled-pair functional) method was used instead. 
The size-extensivity and basis set superposition errors 
(BSSE) were determined to both be around 1-2 cm“^ at 
the n = 7 level of calculation. The calculations were then 
extrapolated to the CBS limit. Despite the meticulous 
care taken in producing this ab initio potential, the re¬ 
sulting equilibrium bond length deviated by 0.13 pm from 
the spectroscopic value. While in excellent agreement for 
a typical ab initio calculation, for a benchmark molecule 
it is a little disappointing that such a powerful calcu¬ 
lation could not achieve a better agreement. Tremen¬ 
dous pains were then mad^^^ to include a host of minor 
corrections to the potential including electron correla¬ 
tion beyond the MR-ACPF method (to a level equivalent 
to full Cl), scalar relativistic effects and diagonal Born- 
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Oppenheimer corrections (DBOC). However, the energy 
effect of each of these additional contributions (all cal¬ 
culated using much smaller basis sets than used in the 
MRCI calculation) was small and mutually cancelling, 
while the improvement in bond length was just 0.04 pm. 
Furthermore, the bond dissociation energy was 90 cm“^ 
higher than the best experimental dataP^ at that point 
(17590 ± 200 cm“^), though smaller than the raw ah ini¬ 
tio result (17699 cm“^). The second papeiP^ concerned 
taking that experimental work by Le Roy et aP^ and fit¬ 
ting a new analytical function to the spectroscopic data. 
Just 5% of the rotationless adiabatic potential has yet to 
be covered experimentally and still the quantum number 
of the highest bound vibrational state is unclear. Again 
great attention to detail was conducted in this study, in¬ 
cluding mass corrected Ce, Cg and Cio constants for the 
different hydrogen isotopes. Despite this, the question of 
whether u = 12 is the highest vibrational level is still un¬ 
resolved but the experimental Dg is now within 1 cm“^ 
of the corrected ab initio value. If confirmed by future 
experiments, this would be an amazing achievement for 
ab initio quantum chemistry. 

Given the difficulties in describing the relatively sim¬ 
ple BeH, the task of completing accurate potentials for 
the heavier system BaH, that contains an additional 52 
electrons, looks formidable. Furthermore, relativistic ef¬ 
fects that were tiny in BeH are amplihed in such a heavy 
hydride. However, 46 or even 54 of the electrons can be 
replaced with an ECP and many of the scalar relativistic 
effects can be included in that potential. The larger core 
was used by Lan^ and earlier by Preuss and co- worker^ 
but in both instances the predicted ground state bond 
length was calculated to be much too short using a triple 
zeta basis set. While Preuss and co-workermanaged 
to achieve 1 pm accuracy using a smaller (4s,4p,2d) basis 
set, such behavior is symptomatic of simple error correc¬ 
tion. Furthermore, both De and oje were much lowei!^ 
than experimental estimates at that time. It should 
not be forgotten that, despite the short bond length, 
the 2-electron, CPP -I- ECP calculation of FC factors 
is sufficiently accurate to be useful for screening good 
laser cooling transitions. For accurate calculation of the 
molecular constants, however, the evidence points to a 
minimum ten-electron basis set for barium. 

Adopting the smaller 46-electron ECP allows the intro¬ 
duction of core-valence interactions with the 5s and 5p 
electrons as well as providing an opportunity to adjust 
the number of electrons used in calculating static and dy¬ 
namic electron correlation. The chaMe in ECP increased 
the predicted re value of Kaupp et aPby a further 10 pm 
but this value was now 4 pm too high. Going one step 
further, Allouche et aP introduced an ^-dependent CPP 
into a Cl calculation of the lowest electronic states of 
BaH using a triple zeta basis set. Although there were 
eleven valence electrons, only three were active in the 
Cl step. The state equilibrium bond length was 

still 3 pm too large, but the agreement is better than 
Kaupp et al^ In addition the molecular constants for all 


the excited states correlating the lowest hve electronic 
states of barium were calculated. They also investigated 
the effect of including spin-orbit coupling in the compu¬ 
tation and demonstrated it made no difference (at the 
100 fm level) to the ground state bond length and al¬ 
tered excited states by < 0.3 pm. However, the spin- 
orbit model adopted was fairly basic and consequently 
the spin-orbit constants were up to 60% too large. The 
calculated dissociation energy, De, of 2.04 eV was signif¬ 
icantly larger than the earlier value of 1.79 eV recom¬ 
mended by FuentealbcP et al from their pseudo-potential 
plus single and double configuration interactive (CISD) 
model, with an equilibrium bond distance of 2.3707 A. 
We note that Fuentealba et aP also obtained a value of 
2.1 eV for De using a pseudo-potential plus local spin 
density (LSD) approximation with an equilibrium bond 
distance of 2.4829 A. That dissociation energy is in excel¬ 
lent agreement with a recent relativistic coupled-cluster 
(RCCSDT) value of 2.062 eV by Skripnokov et aP. This 
is interesting in itself as such agreement between relativis¬ 
tic and essentially non-relativistic calculations is rather 
unexpected in such a heavy hydride. Despite the close 
agreement between the highlighted theoretical bond ener¬ 
gies, they are all larger than the acceptecPi^ experimental 
value. 

Following its discovery in the visible emission spectraP^ 
of the Group H hydrides by Eagle in 1909 the E^H —>• 
0 0 band in BaH was assigned by SchaafsmaP^ 
around the time that Watson began a series of spectro¬ 
scopic studies on BaH recording^ the (I - I) and (2 - 2) 
E^H —>■ X^E+ emission bands (Funke lateil^S added the 
non-diagonal (0 I) and (2 - 1) bands). Watson followed 
that work with the first spectra of the B^E+ —>■ X^E+ 
transitiorP^ in the near infra-red, which revealed an ex¬ 
ceptionally large value for the spin-rotation constant in 
the upper state. As this constant was determined to 
be negative, it suggested the presence of the lower lying 
A^H state that was formally identified by the groujpZESI 
two years later. Also conjectured in these papers was 
the presenc^^ of an even lower lying excited state, 
though this H^A potential was finally directly observecP^ 
by Fabre et al. only in 1987. In 1936 GrunstrbnP^ re¬ 
ported a new higher lying electronic state and analysed 
the (0 0) and (0 - 1) bands of this violet C^E+ —>■ X^E+ 
transition, a similar band having been observed in lighter 
Group H hydrides. Further investigatiorPJ^ revealed that 
the v' = 1 level is abruptly predissociated above iV = 10 
and it was speculated there was a lower bound state with 
an extended equilibrium bond length. 

A series of papers in the 1960 s by Kopp and colleagues 
revealed absorptions in the U\^, identified the D^E+ 
statd^ perturbing the C^E+ vibronic structure, recorded 
the B^E+ ^ X^E+ transition in BaEP^and added addi¬ 
tional vibrational levels to the A^H ^ X^E+ absorption 
spectrgP^ for both hydrogen isotopes. However, due to 
the highly diagonal transitions in the A^H, B^E+ and 
C^E+ ^ X^E+ absorption bands, direct observations of 
vibrational levels were limited to u = 0 - 2 for all the 
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electronic states involved. Only the D^S+ state had a 
potential minimum significantly displaced from that of 
the ground state but has only been studied in absorp¬ 
tion so while the state has been s tudi ed up to v 

= 9, the ground state remains unexplorecP^SI above v = 
3. Furthermore, the often quoted experimental value for 
the dissociation energy in the ground state, Dg < 16350 
cm“^, is based^on the analysis of predissociation in the 
state and assumptions made about the nature of 
mechanism involved. 


II. ^5/Af/r/0 CALCULATIONS 

A. Computational methods 

The initial results reported here were obtained us¬ 
ing a parallel version of the MOLPRCC^ suite of 
ab initio quantum chemistry codes (release MOLPRO 
2010.1). MRCl calculations based on State-Averaged 
Multi-Configuration-Self-Consistent-Field (SA-MCSCF) 
wavefunctionP*^^*^ were conducted to determine the po¬ 
tential energy curves as a function of internuclear dis¬ 
tance r out to bond separation of r = 40 bohr (ao). 
When dealing with dissociating systems, size-consistency 
in the calculation is a particular concern, particularly 
with Cl method^^ restricted to single and double ex¬ 
citations such as the MRCl code in MOLPRO. To es¬ 
timate the missing higher excitations, the a posteri¬ 
ori Davidson correctioipS was applied to all the results 
(MRCl -I- Q). To model the molecular orbitals a variety 
of correlation-consistent Gaussian basis set^^ were tri¬ 
alled; aug-cc-pVnZ, cc-pCVnZ-PP, cc-pwCVnZ-PP and 
aug-cc-pCVnZ-PP. In each calculation, the 46 core elec¬ 
trons of the barium atom were described with an effective 
core potential of the form ECP46MDF from the Stuttgart 
groupl^ leaving the outer ten 5s5p6s electrons. This al¬ 
lows one to test how each basis set performs for the 
molecular constants, e.g.; the equilibrium bond length 
and the dissociation energy Dg while keeping the active 
space constant. 

All the electronic structure computations are per¬ 
formed in the C 2 v Abelian point group symmetry. In the 
C 2 v point group, molecular orbitals are labelled accord¬ 
ing to symmetries ai,bi,b 2 and 02 ; when the molecular 
symmetry is reduced from C^ov to C 2 ,,, the correlating 
relationships are cr —)■ ai, tt —>■ ( 61 , 62 ) , —>■ (oi, 02 ). 
The preferred active space consists of 11 electrons and 
17 molecular orbitals 7ai, 56i, 362 , 202 ; i.e. a (7,5,3,2) 
model. However, to describe the Rydberg character of 
the highly excited states of Ai symmetry it was neces¬ 
sary sometimes to expand the active space to (9ai, 56i, 
362 , 202 ). The averaging process was carried out on the 
lowest five ^ 2 “*", ^H and molecular states of this hy¬ 
dride. 


B. Ab initio potentials using an AV6Z basis 

Typically, the basis set of choice for an MRCl -I- Q cal¬ 
culation including excited electronic states is the aug¬ 
mented correlation consistent type such as aug-cc-pVnZ 
to provide the diffuse Gaussians required to model the 
higher energy electron orbitals. In addition, the larger 
the value of n the more faithful the model wavefunctions 
match the real ones. Therefore, the first computations on 
this hydride were performed using an aug-cc-pV 6 Z basis 
on both atoms (Fig. [^, extending the previous work of 
Lan^. The calculated potentials can be grouped accord¬ 
ing to three broad types: 

(1) Bound potentials based on the interaction between a 

Ba"*" cation and H“ anion e.g. the and A^H states, 

(2) Rydberg states with a single electron held within the 

potential of a BaH+ molecular ion composed of Ba^"*" and 
H“ ions e.g. and 

(3) extended bond length states formed by the avoided 
crossing of largely repulsive valence curves and attractive 
ionic potentials e.g. 

Only doublet states were investigated as the lowest quar¬ 
tet stated are all repulsive. The historical spectroscopic 
assignments were used for the electronic states in Fig. 
[^although it is clear that the adiabatic 3^E+ potential 
is reported as two separate states experimentally (G^E'*' 
and D^E+). The calculated energy gaps at the majority 
of the avoided crossings in BaH are tiny (the case of the 
2 ^E+-X^E+ adiabats being a noteworthy exception) and 
consistent with experiments where, for example, the Ry¬ 
dberg G^E+ vibrational progressiorP^^ is maintained for 
energies well above that of the G/D^E+ states’ closest 
approach (inset. Fig. [^. These AV 6 Z ASE potentials 
are consistent with previous work by Allouche^ et al. 

For the ground X^E+ state the calculated value of Cg 
= 2.18 A is in close agreement to the value of re = 2.16 A 
obtainecP using a barium triple zeta basis set (cc-pV5Z 
basis on the hydrogen), a 56-electron EGP and an i in¬ 
dependent core polarisation potential (GPP). However, 
the accepted experimental valu^ for re is 2.2319 A. The 
AV 6 Z basis set also predicts a dissociation energy Dg 
of 2.2687 eV (18293.08 cm~A while the experimental 
value reported by Kopp et is considerably smaller, 
< 2.0271 eV (16350 cm“^). Both the discrepancies in 
the equilibrium bond length and the dissociation energy 
are rather large for an MRGI-I-Q calculation with such 
a substantial barium basis set. Therefore, the investiga¬ 
tion was widened to include a variety of other basis sets 
in an attempt to improve the agreement between the ab 
initio results for the ground state for this hydride and 
the spectroscopic constants. 


C. Augmented basis sets 

An even tempered augmentation can be used to extend 
both the cc-pCVQZ-PP and the cc-pCV5Z-PP basiP^ 
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bond length r (A) 


FIG. 1. MRCI+Q results for the lowest 5 states of each symmetry ^11 and ^A) as a function of bond length. The 

calculations have been performed with an aug-cc-pV6Z basis set and an ECP with the MOLPRO suite. The potentials shown 
were used to extract values for both the equilibrium bond distance and dissociation energy De. For the ground state, 

the calculated equilibrium bond length is re=2.18 A while the dissociation energy is De = 2.2687 eV. The electronic states of 
the Ba atom fragments are also labelled (only H atoms in the lowest l^S state are produced at these energies). Inset: details 
of the avoided crossing between C/D^E^ states. 


sets. Each basis set is augmented with a single diffuse 
Gaussian function, the augmented exponent a' given by 
the expression. 


a’{£) 


m 


( 1 ) 


where a(£) and /3(£) are respectively the last and second 
last exponent of the appropriate basis set. Calculations 
are performed in the MRCI+ Q approximation for the 
even tempered augmented sets (the additional functions 
are presented in Table M and explicitly optimized aug¬ 
mented cc-pCVnZ-PP ^ = 4 and 5) setd^ (these ex¬ 
plicitly optimized functions are presented in Table [lT|) . A 
large augmented basis set is chosen because in quantum 
chemistry calculations they are known to recover typi¬ 
cally ^ 99 % or more of the electron correlation energy. 


D. CCSD(T) versus MRCI-I-Q 

In the search for an ideal basis set, numerous calcula¬ 
tions with the MOLPRO suite were performed for the 
ground state of this hydride. Table [ITT] presents results 
for the different basis sets used in the MRCI-I-Q approx¬ 
imation, both for the equilibrium bond distance re and 
dissociation energy De, and compares them to CCSD(T) 
calculation^ on BaH. CCSD(T) is currently regarded 
as the most reliable theoretical method for information 
on ground state potentials. In the MRCI-I-Q calcula¬ 
tions, the CASSCF step is performed over the lowest 
three ^Ai, ^Bi, and ^^2 states while only the lowest 
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TABLE I. Augmented functions to the cc-pCVnZ-PP basis 
sel!^ for barium. The even tempered augmentation procee- 
dure is used, where a{£) and P{£) are respectively the last and 
second last exponent in appropriate basis set. 


cc-pCVQZ-PP 

Barium 


Angular momentum function 


Exponent“ 

s 


0.32620 

p 


0.10943 

d 


0.18940 

f 


0.43550 

9 


0.47145 

cc-pCV5Z-PP 

Barium 


Angular momentum function 


Exponent“ 

s 


0.18570 

p 


0.11590 

d 


0.29510 

f 


0.38060 

9 


0.46480 

h 


0.62250 


“Estimated using a'{£) = 


Energy relative to -25.8921314 Hartrees 



FIG. 2. BaH minimum energy in X state, the approach 
to the CBS limit using an active 3-electron model in the 
MRCI-I-Q approximation. Even tempered aug-cc-pCVnZ-PP 
basis sets are used in the calculations. 


TABLE 11. Augmented functions to the cc-pCVnZ-PP basis 
set for barium. For the augmentation proceedure the expo¬ 
nents are optimisecP^ on the energy of the ground state of the 
BaH molecule. 


cc-pCVQZ-PP Barium 


Angular momentum function Exponent 

s 0.0073 

p 0.0060 

d 0.0114 

/ 0.0373 

g 0.0492 


cc-pCV5Z-PP Barium 


Angular momentum function Exponent 

s 0.0071 

p 0.0058 

d 0.0089 

/ 0.0326 

g 0.0393 

h 0.0515 


state appears in the MRCI code. It appears that 
the even tempered augmented basis set aug-cc-pCV5Z- 
PP performs best at MRCI-I-Q level for the ground state, 
matching within 0.1pm the experimental value. The 
explicitly optimized augmented basis set aug-pCV5Z-PP 
gives comparable equilibrium bond distances. In contrast 
to calculations using CCSD(T), the cc-pCVnZ-PP basis 
sets pull the equilibrium bond lengths to much shorter 
values when using the MRCI -(- Q method. Augmenting 
the basis set with diffu se orbitals helps to counteract this 
the CCSD(T), all electron, worls® 


effect. From Table 


III 


provides results (with either the cc-pCV5Z-PP or cc- 
pwCV5Z-PP basis) comparable to the present MRCI-I-Q 
results with the even tempered augmented basis set. The 
MRCI -I- Q calculations are also in excellent agreement 
with the previous RCCSDTI^ and MRCI-I-Q -I- including 
spin-orbital coupling (SOC)pl values for the dissociation 
energy. In addition, this optimal basis set can also be 
used to carry out calculations of excited states, such as 
A^n, a significant advantage of the MRCI method over 
CCSD(T). 

E. Complete Basis Set Limit 

To remove the uncertainty that any close agreement with 
experiment is the result of computational artefacts, it is 
best to extrapolate the energies E„ calculated with the 
VnZ basis sets to the limit of an infinite number of Gaus¬ 
sian functions. For the various basis sets shown in Table 
[nT| the potential energy of the ground state via the com¬ 
plete basis set limit (CBS) is found using the formulae 
derived by JenserP^ with modifications by Karton and 
MartirP^ 


Ecbs = En + 


{Eji — En-l) 

;^ea;p[9(Vn - y/n - 1) - i] 


( 2 ) 


where here n = 5. This expression is used to determine 
all the CBS potential energy curves in this study. Fig[^ 
shows the orderly approach to the CBS energy limit for 
the even tempered augmentation of the cc-pCVnZ-PP 
basiP^ set. While this extrapolation is usually adopted 
for static-correlation only, the small dynamic correlation 
found in the BaH calculations prompted its use here. The 
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TABLE III. Spectroscopic constants for the ground state of the BaH molecule: the equilibrium bond distance re(A) and 

the dissociation energies De in eV and cm“^. Comparison of the present MRCI + Q results with a variety of other reported 
values using correlation - consistent basis set^^ (cc-pVnZ-PP, aug-cc-pVnZ-PP, cc-pCVnZ-PP, cc-pwCVnZ-PP, and aug-cc- 
pCVnZ-PP), both with Cl and CCSDT methods (fc, is frozen core, all, is all-electron) and other previous theoretical work. 
The best experimental values are also included. 


Basis 

Method 

re (A) 

De (eV) 

De (cm ^) 

cc-pVQZ-PP 

MRCI-bQ 

2.2039 

2.2229 

17929.31 

cc-pV5Z-PP 

MRCI-bQ 

2.2329 

2.5961 

20939.44 

CBS 

MRCI-bQ 

2.2376 

2.6580 

21438.71 

aug-cc-pV6Z 

MRCI-bQ 

2.1800 

2.2687 

18293.08 

cc-pCVQZ-PP 

MRCI-bQ 

2.2367 

2.0543 

16569.43 

cc-pCV5Z-PP 

MRCI-bQ 

2.2151 

2.5829 

20832.97 

CBS 

MRCI-bQ 

2.2118 

2.6715 

21547.60 

cc-pwCVQZ-PP 

MRCI-bQ 

2.2367 

2.0587 

16604.92 

cc-pwCV5Z-PP 

MRCI-bQ 

2.2118 

2.5954 

20933.80 

CBS 

MRCI-bQ 

2.2080 

2.6854 

21659.71 

aug-pCVQZ-PP 

MRCI+Q“ 

2.2354 

2.0814 

16787.75 

aug-pCV5Z-PP 

MRCI+Q“ 

2.2332 

2.0871 

16833.78 

CBS 

MRCI+Q“ 

2.2322 

2.0946 

16895.12 

aug-pCVQZ-PP 

MRCI-bQ*’ 

2.2273 

2.1157 

17064.66 

aug-pCV5Z-PP 

MRCI+Q*’ 

2.2204 

2.2384 

18054.33 

CBS 

MRCI+Q*’ 

2.2187 

2.2278 

17968.83 

cc-pVQZ-PP 

CCSD(T), fc" 

2.3302 

2.0311 

16382.30 

cc-pV5Z-PP 

CCSD(T), fc" 

2.3304 

2.0563 

16585.56 

cc-pCVQZ-PP 

CCSD(T), fc" 

2.3266 

2.0494 

16529.91 

cc-pCV5Z-PP 

CCSD(T), fc" 

2.3280 

2.0632 

16641.21 

cc-pwCVQZ-PP 

CCSD(T), fc" 

2.3279 

2.0498 

16533.13 

cc-pwCV5Z-PP 

CCSD(T), fc" 

2.3282 

2.0632 

16641.21 

cc-pCVQZ-PP 

CCSD(T), all" 

2.2378 

2.0298 

16371.82 

cc-pCV5Z-PP 

CCSD(T), all" 

2.2316 

2.0723 

16714.61 

cc-pwCVQZ-PP 

CCSD(T), all" 

2.2379 

2.0355 

16417.79 

cc-pwCV5Z-PP 

CCSD(T), all" 

2.2327 

2.0736 

16725.10 

aug-cc-pVQZ 

RCCSD(T)‘* 

2.2400 

2.0886 

16846.08 

aug-cc-pV5Z 

MRCI+Q-bSOC" 

2.2417 

2.0860 

16825.11 

Triple-zeta 

CIPSI (ECP-bCPP)-'^ 

2.2620 

1.9700 

15889.49 


Pseudo-potential -|- LSD® 

2.4280 

2.1000 

16938.03 


Pseudo-potential -|- CISD® 

2.3707 

1.7900 

14437.66 

Experiment 


2.2319'“ 

<2.0271' 

<16350® 

cc-pCVQZ-PP 

MRCI-bQ (3e)'“ 

2.23670 

2.0543 

16569.43 

cc-pCVQZ-PP 

MRCI+Q (5e)' 

2.23725 

2.1359 

17227.59 


“MRCI-fQ, even tempered basis set augmentation functions, present work. 

*’MRCI-|-Q, explicitly energy optimized augmentation functions, present work. 

'^CCSD(T), fc and CCSD(T), all : private communication from Peterson (2015). 

'^RCCSD(T), Skripnikov and co-workerP. 

"MRCI+Q -b SOC, Gao et aP. 

■'’CIPSI Allouch(P. 

®Pseudo-potential -|- local spin density (LSD)|^, Pseudo-potential -|- CISEP. 

^Experiment, Ram and BernatiP. “Experiment, Kopp et aP. ■’Assuming Ba(®D 3 ) limit. 

*’MRCI-|-Q, active 3-electron model, present work. ^MRCI-|-Q, active 5-electron model, present work. 
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close agreement with the experimental le value found in 
the CBS potential supports this choice. 

In Table |III| the results for the equilibrium bond 
distance rg and dissociation energy Dg from these 
CBS limit, MRCI+Q calculations are compar ed w ith 
the coupled cluster appromixation CCSD(Tj3nillHl2l^ 
in its various imple mentations and with other theo¬ 
retical methodj ^ l ^ l ^^l The density functional theory 
(DFT) work of Wang and Andrew^ using a 6-311-1—|-G 
(3d/, 3pd) basis yielded a value of 2.2520 A and second or¬ 
der Mpller Plesset perturbation theory gave 2.2440 
respectively for the equilibrium bond distance. Gao et 
have conducted calculations for the ground X^I]+ 
and the first excited state A^II of this hydride at the 
MRCI-fQ-fSOC level. Their ab initio results for the 
spectroscopic constants (2.2417 A) and Dg (2.0860 eV) 
of BaH (X^E+) ground state are similar to those ob¬ 
tained from the relativistic - coupled-cluster result^ 
(RCCSDT), respectively 2.24 A and 2.0840 eV, and in 
reasonable agreement with the available experimental 
result^^. The other ab initio results of great accuracy are 
those using CCSD(T) where the cc-pCV5Z-PP (all elec¬ 
tron) result is within 0.02 pm of the experimental valu^ 
but if extrapolated to the CBS limit this agreement 
will be somewhat poorer. Using an aug-cc-pCVnZ/CBS 
MRCI-I-Q calculation the ab initio equilibrium bond dis¬ 
tance here is within 0.031 pm of the best spectroscopic 
value. 

III. MODELLING THE GROUND STATE 

A. Spectroscopic constants of the ground state 

An important requirement is ensuring that the final po¬ 
tential behaves correctly at all interatomic distances, 
particular at long-range where dispersion forces are 
dominanlP^. In atomic units, the asymptotic atom-atom 
potential has the forrJ^ 

= + (3) 

where Ce, Cg etc. are the polarizabilities with Ityo the 
asymptotic limit (atomic products), r is the atom-atom 
internuclear separation. The dynamic polarization of the 
ground state, which allows the computation of the lead¬ 
ing Ce term (148 a.u.), has been calculatecP^ by Dere- 
vianko et al. To calculate the following Cg term, the 
formula derived by Tan^^ requires the quadrupole po¬ 
larizabilities but there are two theoretical valuei^^*^ for 
Ba that produce Cg values that differ by more than 60%. 
The smaller value, 11242 a.im is consistent with the cor¬ 
responding term calculatecP^ for the neighbouring hy¬ 
dride CsH (11710 a.u.), and is therefore the value recom¬ 
mended at present. There is no theoretical or experimen¬ 
tal value at present for Cio in BaH so the expansion was 
truncated at Cg: the corresponding calculated valu^^in 
CsH is 1.194 X 10® a.u. 



FIG. 3. Analysis of ab initio ground state of BaH using 

LEVEL^. The ab initio data was taken with the aug-cc-pCVnZ 
basis sets (n = 4 and 5) and extrapoliated to the CBS limit. 
The vibrational levels fonnd in the fit are also presented. 



FIG. 4. Comparison of the gronnd state energy for the BaH 
molecnle in the region of the equilibrium bond distance us¬ 
ing both an active 3-electron and 5-electron model in the 
MRCI-I-Q approximation. The cc-pCVQ-PP basis seP^ is 
used for both calculations. The absolute energies are adjusted 
by 25 E/r. 

The LEVEL^^ program (Version 8.2) was used to cal¬ 
culate the rovibrational levels of the resulting ah initio 
X^E+ aug-cc-pCVnZ/CBS potential. A smooth poten¬ 
tial was produced through interpolation of the ab ini¬ 
tio points with a 4-point piecewise polynomial (NUSE=4, 
IR2=1), and extrapolation past r = 15.88 A using a 

fixed Cq = 7.110512 xlO® cm“^ A ® (no Cg or Cio) 
(ILR=2). Based on this potential the equilibrium bond 
length and dissociation energy were determined to be 



























TABLE IV. Vibrational term energies in the ground state of the BaH molecule. Comparison of the present MRCI+Q 

results with the aug-cc-pCVnZ-PP even tempered basis set taken to the CBS limit and experimental data. All energies in 
cm“^. 


v" 

MRCI 

Relative value“ 

ExperimentaP 

RKR*" 

AE RKR 

0 

581.3491 

0 

0 

580.53 

0.8191 

1 

1722.2738 

1140.9247 

1139.2896 

1719.84 

2.4338 

2 

2832.9237 

2251.5746 

2249.6062 

2830.15 

2.7737 

3 

3915.9012 

3334.5521 

3331.1192 

3911.46 

4.4412 

4 

4981.7132 

4400.3641 


4963.77 

17.9432 

5 

6024.9686 

5443.6195 


5987.08 

37.8886 

6 

7029.1845 

6447.8354 


6981.39 

47.7945 

7 

8006.0703 

7424.7212 


7946.70 

59.3703 

8 

8957.6993 

8376.3502 


8883.01 

74.6893 

9 

9878.6102 

9297.2611 


9790.32 

88.2902 

10 

10772.2217 

10190.8726 


10668.63 

103.5917 

11 

11633.1528 

11051.8037 




12 

12461.4104 

11880.0613 




13 

13257.3637 

12676.0146 




14 

14016.6614 

13435.3123 




15 

14729.2373 

14147.8882 




16 

15386.3161 

14804.9670 




17 

15966.4205 

15385.0714 




18 

16436.4095 

15855.0604 




19 

16749.4237 

16168.0746 




20 

16882.9819 

16301.6328 




21 

16894.9472 

16313.5981 





“E(MRCI) - E(experimental) = AE (u' - v”) cm (1 - 0) = +1.64; (2 - 1) = +0.33; (3 - 2) = +1.46. 
*’Based on experimental data prioi^to 1982. 



FIG. 5. Tuning r^e/ to improve the fittiM of the ab initio 
ground state of BaH using betaFn^^. 


2.2322 A and 16895 cm“^. Twenty-two vibrational levels 
were predicted for the system, with z; = 21 lying just 0.2 
cm“^ below the dissociation limit. Comparison of the 
relative energies of the first four vibrational levels based 


on this potential with experimental val uej ^ show good 
agreement, all within 1.7 cm“^ (Table IV). There is a 
significant jump in the error with respect to the earlier 
RKR levels at u = 4, a result of the lack of experimental 
data for the higher vibrational levels and consequently 
the present results are probably more reliable for those 
missing vibrational levels. 


B. 3 electron versus 5 electron MRCI+Q 

Though the calculated bond length is only <0.031 pm 
away from the experimental value, another effort was 
made to explore possible improvements to the accuracy 
of the ab initio potential. In order to quantify the ef¬ 
fect of including additional dynamic electron correlation, 
further MRCI+Q calculations were performed with the 
cc-pCVQZ-PP basis seiP^ using both 3- and 5-electrons 
active at the Cl stage. From the ab initio work shown in 
Fig. ID it is found that the active 5-electron model lowers 
the absolute (total) energy of the ground state minimum 
for this hydride by just 3 milli-Hartrees (mE/j). The 
subsequent change in the equilibrium bond distance Ve 
was also tiny, an increase of 0.055 pm. These very small 
changes in the total energy and equilibrium distance sub¬ 
sequently justifies the use of the 3-electron model for 
both the ground and excited states considered as opposed 
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to the more computationally intensive active 5-electron 
model. 


C. Including spectroscopic data 


There is an extensive selection of pair-potentials in the 
literature to describe the r dependence of the interaction 
energy in a diatomic molecule (see the supplementary in¬ 
formation from Xie et al for a comprehensive listP^ with 
over a hundred examples). Two of the most success¬ 
ful functions are the Tang-Toenniej^ potential (and its 
variants) and the Morse-Long-Range^ (MLR) function. 
The ab initio data here was fitted to a MLR potential 
using betaFIT (version 2.1)1^ in order to find an analyt¬ 
ical function to express the potential. Prior to fitting, 
the potential was shifted by 0.03 pm so that the ab initio 
potential minimum was aligned with the experimental re 
value to an accuracy better than 0.01 pm. This simplifies 
the inclusion of experimental data at a later stage. The 
MLR potential is described by 

RMLR(r) = Se fl - (4) 

V u[re) J 


where 2)e is the dissociation energy and u{rf.) is value of 
u{r) at the equilibrium bond distance. 

Both Vmlr and /3(r) are expanded in terms of two 
internuclear distances rg and r^e/, according to the radial 
variable 


Vnir) = 


-Fr” 


( 5 ) 


where r^ = or r^ef and the latter is usually greater 
than or equal to r^. As the fit is viable for a large range 
of Vrefi it is worth tuning this value to improve the fit 
(Fig. MLR functions producing local minima were se¬ 
lected as starting points for later fitting to spectral data. 
Typically, two different values of n are required for an 
acceptable fit (in the present work, n = 8 and 3) in the 
/3(r) expansion. The first issue is what to set as the error 
in the ab initio points. Since the vibrational frequencies 
are calculated to within an average of 2 cm“^ the error in 
each ab initio point was initially set rather arbitrarily as 
10 cm“^ (5 X 10“^ E?i). However, Dattani and Le Roj^^ 
have shown that even in the 6e Li 2 molecule a deviation 
in De of up to 68 cm“^ still persists in the calculated 
b^n state. Therefore, an increased error of 20 cm“^ was 
adopted as a compromise value. The final analytical po¬ 
tential is reported as MLRp(q^ {Np), where p and q are 
the two values of n required and consists of Np expansion 
coefficients. A similaJ^ potential (replacing an earlier 
versioiP^ that used only one expansion co-efficient p) has 
recently been used to describe the ground state of MgH, 
the only hydride other than H 2 with an experimental 
dissociation energy determined to sub-wavenumber ac¬ 
curacy. The fitted Se is just over 4 cm“^ lowe r than the 
result from the ab initio potential (Table [VT| ), as it was 
allowed to vary slightly during the fitting process. 


The parameters fitted to MLRg g"^ (12) from the ab 
initio potential were then transferred to 
(version 2 .0) and combined with B^I]+ —)■ 

emissioiP^l2lIlll x^E+ infrared dateP^. The earlier 
emission data was quoted to 0.01 cm“^ while the stated 
resolution of the infrared data was 0.005 cm“^. 

While Tref = 3.34 A does not correspond to a local 
minimum in Fig. those MLRs obtained from the local 
minima in betaFIT were not found to be as successful. 
p = 8, q = 3 were selected as the best parameters af¬ 
ter screening other combinations. As DPotFit utilizes 
a dimensionless root-mean-square-deviation (dd) as its 
means of determining how closely a model matches the 
observed data, it is crucial that a reasonable assessment 
of the uncertainty in the observable data is made. 


dd = 


-^data 

i ^ 

^ -^calc(^) -^obs(^) \ 

A^data 

V ^obs(^) / 


( 6 ) 


Thus, as the uncertainty in the ab initio Fobs was set to 
20 cm“^, deviations from the fitted MLR do not lead to 
a large increase in the dd. At the chosen Tref the dd be¬ 
tween the MLRg-3^(12) potential fitted by betaFIT and 
the ab initio points was only 0.4570 representing a close 
fit within the selected uncertainty. By contrast the origi¬ 
nal emission data was quoted to much higher accuracy, so 
any deviations between the model and the spectral mea¬ 
surements can lead to a significant increase (left hand 
data. Table VI) in dd. Due to the arbitrary nature of the 


uncertainty assignment in the ab initio data, the associ¬ 
ated dd should not be compared directly to the dd in fits 
to experimental data. 

As DPotFit tries to fit to very accurate spectroscopic 
data, the initial uncertainty in the fit is very high. Al¬ 
lowing j3 parameters to become fitted parameters pro¬ 
duces a final potential which matches the experimental 
data much more closely than the ab initio potential alone. 
During the fitting, both Cq and Cg were fixed at the rec¬ 
ommended theoretical values, and were incorporated as 
the long range tail of the MLR potential with Douketis- 
type dampin^^with s = — 1. Additionally, the effect of 
spin rotation coupling is incorporated through a modifi¬ 
cation to the MLR potential^ using the expansion coef¬ 
ficients listed in Table lYul 

As the considered spectral data covers a range of 
isotopomers of BaH (electronic transitions of ^^®BaH 
and ^®®BaD; infrared transitions of ^®®BaH, ^®^BaH, 
^®®BaH, ^®®BaH), all of which have very small val¬ 
ues of reduced mass, it is valuable to consider the 
impact of the Born-Oppenheimer breakdown (BOB)l^ 
corrections. DPotFit incorporates BOB corrections 
through the atom-dependant potential, with the effects 
parametrised within the radial strength functions: 




Stiir) = vZ (•■) + [1 - y?Jr)] X < «z (0* (7) 


2=0 















10 


TABLE V. Summary of the spectral data used when fitting with DPotFit. A total of 1769 B^E-X^E transition and 409 X^E'*" 
IR emission lines were incorporated. 


Year 

Transition 

Isotope 

V V (t/max) 

Lines 

Unc / cm ^ 

193#^ 

B2e+ X2E+ 


0-0 (37.5) 

112 

0.01 

193p^ 

B2E+ ^ X2E+ 


1-0 (41.5) 

142 

0.01 




1-1 (31.5) 

106 

0.01 




2-1 (37.5) 

126 

0.01 

196d^ 

B2E+ ^ X2E+ 

^®®Ba2H 

0-0 (52.5) 

197 

0.01 




1-0 (47.5) 

177 

0.01 




0-1 (42.5) 

146 

0.01 




1-1 (50.5) 

179 

0.01 




2-1 (46.5) 

162 

0.01 




1-2 (33.5) 

94 

0.01 




2-2 (44.5) 

142 

0.01 




3-2 (39.5) 

125 

0.01 




3-3 (29.5) 

61 

0.01 

199^2^ 

IR 


1-0 (28.5) 

94 

0.005 




2-1 (20.5) 

80 

0.005 




3-2 (29.5) 

77 

0.005 




1-0 (22.5) 

59 

0.005 




1-0 (22.5) 

59 

0.005 




1-0 (19.5) 

40 

0.005 


TABLE VI. The MLRpy|^(V^) potential parameters obtained throngh a fit to the ab initio potential, and after initial optimi¬ 
sation with spectral data using DPotFit. All energies in cm“^. 


MLR|;3'^(12), ab initio only |MLR|J"^(12), ab initio + spectra 


De 16891.5947 
re 2.2319 A 
dd 266.447 

De 16841.8651 | 

re 2.2318 A 
dd 6.083 

Po 

0.114040 

07 

-87.820760 

Po 

0.096691 

07 

-86.606654 

Pi 

-7.619529 

0s 

-171.874803 

Pi 

-7.616896 

Ps 

-154.070717 

02 

-17.121381 

09 

-5.387603 

P 2 

-16.467554 

P9 

34.3428882 

05 

-19.233860 

Pio 

256.820563 

P3 

-16.902752 

PlO 

324.143391 

Pi 

3.673494 

Pii 

260.900999 

Pi 

6.565879 

pll 

334.502110 

/35 

44.305094 

P 12 

83.391642 

P 3 

44.183256 

P 12 

116.057107 

06 

31.801460 



06 

28.715007 






Rtir) = yllSr)tt + [1 - vZSr)] (8) 


z=0 


Where is the atom-dependant adiabatic BOB ra¬ 

dial strength function and R^Z''") the non-adiabatic BOB 
radial strength function. Incorporating fitted terms (Ta¬ 
ble VII) for iVad = Nna. = 2 along with fitted MLR pa¬ 
rameters (Table VI) while keeping re fixed at 2.2319 A 


produced an MLR potential with dd = 6.143. Allowing 
Te to also be a fitted parameter lowered Sg by nearly 
50 cm“^ but further reduced dd to 6.083. While sub¬ 
stantially higher than what would usually be deemed a 
“good” fii!^ {dd < 1 ), comparison with results from the 
most recent studj^^of the E(^n) —X(^S+) transition, a 
data set not used in the DPotFit analysis, highlights the 


much improved agreement with the experimental obser¬ 
vations (Table |VIII) than the raw ab initio results from 
Table 13 

In addition, a number of further details were evident 
while developing the hybrid potential. The experimental 
value^^for the sp in-rot ation constant 'y{v) follow a very 
clear v dependenc^^SEl (Fig. that can be modelled as 
a simple Taylor (linear) expansion; 


7(^^) = 70 + 7i + 0 (9) 

with 7 o = 0.19456 cm“^ and 71 = —4.9813 x 10“^ cm“^. 
Finally, the experimental dissociation energy quoted in 
the literature has its origin in the sudden onset of rapid 
predissociation in the v' = 1, N' = 10 rovibronic 
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V 



FIG. 6 . The v dependence of the experimental spin-rotation 
constants 7 from the experimental worlP^ of Walker, Hed- 
derich and Bernath. 


TABLE VII. Additional parameters defining the spin- 
rotation coupling and the BOB corrections in the final MLR 
potential, obtained after initial optimisation with spectral 
data using DPotFit. 


MLR gj^(12). Additional parameters 




3 



Wo 

0.058666 



Wi 

-0.034468 



wf 

-0.015660 


Pad 

6 

Vad 

6 

ur 

14.158630 

Mo 

-81.921108 

uf^ 

-106.246129 

M? 

244.705054 

nr 

-3.266579 

M 2 

-25.518624 

wSf 

0.000000 

mS) 

0.000000 

Pna 

3 

Qna 

3 

>Ba 

^0 

-0.007379 

to 

0.007367 

>Ba 

0.003445 

tf 

0.000210 

>Ba 

^2 

-0.015764 

^2 

-0.049889 

>Ba 

^00 

0.000000 

^00 

0.000000 


level. This level lies 25942.6 ± 0.5 cm“^ above the 
minimum. Assuming that this indicates the dissociation 
limit of the (diabatic) potential lies above 

v' = N' = % (25874.6 ± 0.3 cm“^) and identifying the 
asymptote as belonging to Ba(^Di) and not Ba(^D 3 ), this 
would fix the dissociation energy Dg of the ground state 
between 16842.6 cm“^ and 16910.6 cm“^ which is con¬ 
sistent with the present ab initio value and other recent 
theoreticaP results, while the DPotFit result is within 
1 cm“^ of the lower limit. It is therefore tentatively sug¬ 
gested that the correct dissociation limit for the 


FIG. 7. MLRg' 3 '^( 12 ) fit to ab initio data points, and final 
MLRgJ^(12) fit to experimental data. 

state is actually Ba(^Di) -F H(^Si) and at long range 
ought to be labelled Di(l) in Hund’s case (c). 

D. Excited states 

Extending the ab initio study to excited states, the lowest 
energy state was calculated usi^ the even tempered 
augmented cc-pCV5Z-PP basis selP^. Despite being car¬ 
ried out with a CASSCF calculation including 5Ai, 5Bi 
and 5 A 2 states and the MRCI being carried out over 
5Bi states, the resulting potential has an le of 2.2698 A, 
within 0.38 pm of the best experimental value from Bar- 
row (quoted in Allouch^ et al). The calculated differ¬ 
ence in and A^II bond lengths is therefore Are 

= 3.97 pm compared to the experimental difference of 
4.16 pm, a discrepancy of just 2 mA. Of course, for 
a truly faithful reproduction of the A^II state not only 
must spin-orbit coupling be included in the calculation, 
but also the angular momentum coupling with other elec¬ 
tronic states. However, for FC factors the fact that Are 
is so close to experiment is reassuring and small energy 
shifts in general do not have a large effect. 

The two potentials were then used to obtain FC factors 
between the X^E+ and A^H vibronic states using LEVEL. 
The highly diagonal nature of the A^H ^ X^E+ tran¬ 
sition is clear when the wavefunctions from both states 
are compared as in Fig. The diagonal FC factor /oo 
is almost identical to that found in the earlier study by 
Lan^ using the larger ECP and a much smaller triple 
zeta quality basis set, though the agreement reduces with 
each increase in vibrational quantum number (by /as the 
difference is 5.8%, compared to 0.1% for /go). Using 
the revised dissociation energy Dg of the ground state 
(< 16910.6 ± 0.5 cm“^), the separation in bar- 
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TABLE VIII. Vibrational term energies in the ground 

state of the BaH molecule using DPotFit. Comparison of 
the vibrational levels from the present combined ab ini¬ 
tio/spectroscopic potential and experimental^ data indepen¬ 
dent of fit. All energies in cm“^. 


V 

Experiment 
(Ram et at) 

DPotFit“ 
re fix 

AE 
re fix 

DPotFit'’ 
re fit 

AE 
re fit 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 

1139.2896 

1139.2884 

-0.0012 

1139.2897 

0.0001 

2 

2249.6062 

2249.6084 

-0.0022 

2249.6078 

-0.0016 

3 

3331.1192 

3331.1212 

-0.0020 

3331.1234 

-0.0042 


“ Zero Point Energy = 580.6050 cm ^ 
Zero Point Energy = 580.5910 cm~^ 


ium (9372 cm“^), the Term energy of the minimum of 
the A^n potential (excluding spin-orbit separation) from 
Barrow^ and co-workers and finally the X^S+ zero point 
energy, the experimental dissociation energy De of the 
A^n state is < 16001.3 ± 2.0 cm“^ which is in good 
agreement with our calculated value of 16076.75 cm“^. 
This is further support for our proposed revision of the 
experimental value for the BaH dissociation energy 

De to < 16910.6 ± 0.5 cm“^. 

Using the same basis set, ab initio calculations were 
performed in the MRCI-I-Q approximation for five states 
of each of the ^Ai, ‘^Bi and ^^2 symmetries for bond 
length out to 40 Uq. The results from all these calcula¬ 
tions are presented in Fig. The ^Ai states are the 
most common symmetry at lower energies and therefore 
valence states dominate in Fig. Comparing to the 
previous AV6Z calculation, the less diffuse nature of the 
largest Gaussian functions in our adopted aug-pCVnZ- 
PP basis sets has resulted in the greater weighting of 
valence states over Rydbergs and so the state has 

disappeared. To describe the Rydberg character of the 
states, it’s necessary to expand the active space in 
the MRCU-Q calculations to at least (8ai, 5bi, 3b2, 22 ). 


IV. SUMMARY AND CONCLUSION 

In this paper we have attempted to develop a reliable 
potential energy curve for the X^E+ state of the BaH 
radical. Encouraged by the close agreement in the lit¬ 
erature between relativistic and non-relativistic ab ini¬ 
tio calculations, we have used the MRCU-Q method to 
determine the shape of the X^E+ potential and the dis¬ 
sociation energy. With an even-tempered augmented cc- 
pCVnZ-PP basis set taken to the CBS limit, we have cal¬ 
culated the equilibrium bond length to within 0.031 pm 
of the experimental value and of comparable accuracy 



FIG. 8. The almost identical first four vibrational wavefunc- 
tions in the X and A states of BaH. The larger equilibrium 
bond length displaces the excited state slightly but the overlap 
is sufficient to ensure diagonal Franck-Condon factors close to 
1. Also labelled are the diagonal Franck-Condon (FC) factors 



FIG. 9. Low lying states of BaH using an active 3-electron 
model in the MRCI-I-Q approximation. The even tempered 
aug-cc-pCV5Z-PP basis set is used for the calculations. 

to the more trusted CCSD(T) technique. While the ten 
5s5p6s electrons from barium are necessary to calculate 
the static correlation energy, only the valence 6s electrons 
are needed to cover the majority of the dynamic correla¬ 
tion. The computed dissociation energy (16895.12 cm“^) 
is in good agreement with other theoretical values sug¬ 
gesting a revision of the experimental dissociation energy 
to De <16910.6 cm“^. The calculated dissociation en¬ 
ergy of the excited A^H state (16076.75 cm“^) is also 
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consistent with this experimental value. The vibrational 
levels of the ab initio potential achieve an average 

agreement of just over 1 cm“^ with the best experimen¬ 
tal values. Refining the potential with experimental data 
using DPotFit improves the agreement to better than 
0.005 cm“^ for the lowest vibrational levels. 
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